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Abstract. We consider the interaction of a slowly rotating unmagnetized neutron star with a hot (ion supported, 
ADAF) accretion flow. The virialized protons of the ADAF penetrate into the neutron star atmosphere, heating 
a surface layer. Detailed calculations are presented of the equilibrium between heating by the protons, electron 
thermal conduction, bremsstrahlung and multiple Compton scattering in this layer. Its temperature is of the 
order 40-70 keV. Its optical depth increases with the incident proton energy flux, and is of the order unity for 
accretion at 10 _2 -10 _1 of the Eddington rate. At these rates, the X-ray spectrum produced by the layer has a 
hard tail extending to 100 keV, and is similar to the observed spectra of accreting neutron stars in their hard 
states. The steep gradient at the base of the heated layer gives rise to an excess of photons at the soft end of 
the spectrum (compared to a blackbody) through an 'inverse photosphere effect'. The differences with respect 
to previous studies of similar problems are discussed, they are due mostly to a more accurate treatment of the 
proton penetration process and the vertical structure of the heated layer. 
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1. Introduction 

The spherical accretion of matter onto the surface of 
an unmagnetiz ed neutron star (NS) has att racted much 
attention since Zel'dovich fc Shakura (1969 ) [henceforce 



ZS69] first addressed this problem more than thirty years 
ago. As a model for the then newly discovered X-ray stars, 
these authors considered a neutron star whose surface was 
heated by radially infalling gas, and modeled this gas as 
consisting of freely falling ions. A hot X-ray emitting layer 
is formed, the temperature of which depends on the pen- 
etration depth of the protons into the atmosphere of the 
star. The model could reproduce the observed hard X- 
rays, but was eclipsed by accretion disk models once it 
was realized that accretion is not radial because of the 
angular momentum constraint. The disk model, however, 
only explains the sub-keV part of NS spectra, somewhat 
in contradiction with observed spectra. In addition to the 
optically thick disk, an additional source of hot gas thus 
had to be found to produce the hard X-rays. This was 
proposed in the form of the so-called two-temperature ac- 
cretion flows ( Shapiro et al. 1976]; jtchimaru 1977 ; Recs 



et al. 1982| ; |Narayan fc Yi 1994 |1995| ). In this iOTrh of ac 



cretion, now called ADAF, the flow is geometrically thick, 
supported by virialized protons, while the electrons stay 
at a lower temperature around 100 keV due to their high 
radiative energy losses and the low rate at which Coulomb 
interactions transfer accretion energy from the ions to the 



electrons. While it is not entirely clear at the moment how 
such a flow would originate from the cool disk present at 
larger distances from the star's surface, the observed hard 
X-ray component is often taken as an indication of its ex- 
istence. If accretion near the star takes indeed place in the 
form of such a two-temperature plasma, the star's surface 
is exposed to protons with energies around 50 MeV, and 
the physics is very much like the model proposed by ZS69. 
It is therefore interesting to pick up this line of theory and 
revisit proton-injection models in context of NS-ADAF 
accretion. 

In ZS69 the importance of the penetration depth of 
the accreted protons for the outcoming spectrum was rec- 
ognized. The penetration length in terms of the amount of 
material required to stop the protons was estimated and 
used as a model parameter. The injected energy of the in- 
falling protons per unit time and mass in the neutron star 
atmosphere was then uniformly distributed in this stop- 
ping layer. A second parameter of their model was the 
accretion rate (or luminosity). Their computed spectra es- 
sentially show a blackbody spectrum with a high-energy 
tail due to Comptonization in the heated atmosphere. 

A more detailed numerical approach was presented by 



Alme & Wilson (1973). They introduced more physics into 
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their model and also realized that the proton deceleration 
depends on the atmospheric temperature. They solved the 
atmospheric structure and the beam deceleration simulta- 
neously in a time dependent evolution of the model. But 
they did not account for the strong dependence of pro- 
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ton deceleration on the local proton velocity and used an 
estimate only for the amount of material needed to stop 
the infalling protons. Their resulting spectra again show 
a blackbody plus an additional high energy tail. 

Bildsten et al. (1992) were interested in the fate of the 
accreted CNO elements in the neutron star atmosphere. 
They assumed an isothermal atmosphere at a tempera- 
ture defined by the accretion rate and the proton stopping 
depth. Again a similar estimate for the penetration depth 
of the ions was used. 

The model o f spherical accretion onto a neutron star 
was revisited by Turolla et al. (1994 ) . They confirmed the 
existence of the form erly known "cold" solutions (as in 
Ahnc &: Wilson 1973 ). In addition to this, they also found 
"hot" solutions for the same luminosities. But their work 
was limited by the fact that the spectrum was just de- 
scribed in terms of a mean photon energy. 

Zane et al. (1998) confirmed the existence of the hot 
solutions with a much more accurate treatment of the ra- 
diative processes including pair processes. The emergent 
spectra of the hot solution should be characterized by hard 
spectra peaking at ss 100 keV. Both models use the stop- 
ping depth as a free input parameter. The effect of thermal 
conduction was not included in these calculations. 



In continuation of this series of papers Zampieri 



et al. (1995) presented solutions for neutron stars accret 
ing at low rates, i.e. with luminosities Lns — 10~ 7 - 
10~ 3 -LEdd- They showed that the deviations from a black 
body spectrum with the effective temperature of the neu- 
tron star increase with decreasing accretion luminosities. 

Proton illumination is not only restricted to the accre 
tion of matter onto the surface of a neutron star. Spruit 



(1997 ), |Spruit fc Haardt (2000|) and [Dcufcl fc Spruit 



(2000] ) have shown that this process might also be of im- 



portance for accretion disks embedded in a hot corona. 



With an improved version of the code used in Deufel & 
Spruit (20001 ) we reconsider the accretion of matter onto 



the surface of a neutron star (i.e. ion illumination of a 
neutron star surface). The code now includes a better 
treatment of the radiative processes by including ther- 
mal emission due to bremsstrahlung. We also allow for 
energy redistribution due to electron thermal conduction 
within the atmosphere. The interaction between the ac- 
creted protons and the neutron star atmosphere and the 
radiative processes are computed time dependently in a 
one-dimensional, plane-parallel approximation. The den- 
sity distribution through the atmosphere is found from 
hydrostatic equilibrium including the pressure from the 
penetrating ions. 

In the following we give a complete description of our 
model. First the heating of the electrons via Coulomb in- 
teractions is presented. We show that the use of a stopping 
depth treated as a free parameter is not sufficiently accu- 
rate for the calculation of the proton stopping. The radia- 
tive processes are treated by solving the radiative transfer 
equation instead of the Monte Carlo method used before. 
This has advantages in terms of computation speed and 



the ability to treat layers of arbitrary absorption optical 
depth. 

2. The Neutron star illumination model 

2.1. Proton stopping by Coulomb collisions 

We consider protons in the vicinity of a neutron star of 
mass Af* and radius In ADAF accretion the typical 
energy of such protons at the surface of a neutron star is 
the virial temperature, 



E p k T v [ r 



GM*m P 
3i?* 



46- 



MeV /M*\ /10km\ 



proton \ M & J \ i?* J 



,(1) 



where G is the gravitational constant, m p the proton mass 
and k the Boltzmann constant. 

If protons with such energies encounter a much a cooler 
and dense medium, they will be stopped very efficiently by 
Coulomb interactions with the cold electrons inside their 
Debye sphere. The stopping of a fast particle in an ionized 
plasma was quantitatively discussed by Spitzer (1962). In 
a time di, such a proton with energy E p will loose an 
amount of energy dE, 



dE . . 

= -e{E p )n c v p , 



(2) 



where n c is the electron density of the neutron star atmo- 
sphere and Vp w 3kT/m p is the proto n velocity. s(E p ) is 
called the stopping power of a plasma. Ryter et al. (1970| ) 
show that it is given through Spitzcrs's formalism by 



e(E) 



47re 4 



m c v„ 



In A [ip(x) - xip'(x)]. 



(3) 



Here InA = ln[(3/2e 3 )(fc 3 T c 3 /7rn c ) 1 / 2 ] is the Coulomb 
logarithm, <tt is the Thomson cross-section, ip( x ) an d 
ip'{x) are the error function and its derivative, and x = 
(m e v p /2kT e ) 1 / 2 is the ratio of the incident proton velocity 
to the thermal velocity of the electrons. This formula holds 
for non relativistic electron and proton temperatures. 

Now we can express the energy loss of a proton moving 
at an angle 9 with respect to the vertical coordinate z, per 
unit of vertical Thomson depth, dr — n c <jT<iz, by 



dE 
dr^ 



47re 4 



or rn e v p cos f 



\nA[ip(x) - xtp'(x)]. 



(4) 



In a time-dependent evolution of the model we fol- 
low the protons numerically through the neutron star at- 
mosphere with temperature T(t) and density n c (r). We 
record the loss of kinetic energy of the protons as a func- 
tion of optical depth. This yields the local time dependent 
heating rate A p (r). In this way the stopping depth of the 
protons is obtained self consistently. 

As mentioned in Sect. |l| the penetration depth signif- 
icantly influences the outcoming spectrum of the proton 
illumination model. Usually a global estimate is used to 
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Table 1. Estimated stopping depth (est) and calculated 
stopping depth (num) [according to Eq. (|J)] in units of 
Thomson optical depth for a isothermal atmosphere with 
temperature T e and uniform electron density n c = 10 22 
cm~ 3 for various proton energies. 







Temperature T e [keV] 


T p [MeV] 


1 




10 


100 




num 


est 


num est 


num est 


23 


0.2 


0.8 


0.4 0.8 


6.8 7.8 


46 


0.8 


3.1 


0.9 2.4 


10.2 12.4 


92 


3.1 


12.4 


2.6 8.9 


15.9 21.6 



calculate the stopping depth of the protons in an atmo- 
sphere. The estimate for the stopping length A s is obtained 
following Spitzer (1962), 



As — Vr> t s 



(5) 



where t s is the slowing down time scale given by [Spitzer 



(1962))] Eq. (5-28). If we express the stopping length in 
terms of Thomson optical depth, we get 



m c 3 In A ip — xtp' ' 



(6) 



where v z is the vertical component of the proton ther- 
mal velocity and c is the speed of light. This expres- 
sion is inaccurate because of several factors. Most impor- 
tantly the proton velocity enters this formula as the fourth 
power. As the proton velocity decreases, this factor drops 
rather fast and shortens the penetration depth consider- 
ably. Counteracting is the function ip — xip'. For x > 2 this 
function is unity. Bu~E tor U < x < 2 it is small, mcreas- 
ing the stopping lengtli but also reducing tlic energy loss 
[cf. Eq. (^)J. Even the Coulomb logarithm should not be 
taken as constant (as it is often done). In the wide range 
of electron temperatures and densities of interest (as solu- 
tions from e.g. |Turolla et al. (1994 ) suggest) the Coulomb 
logarithm can easily vary by a factor of two. 
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Fig. 1. Energy loss rates of a hot thermal distribution of 
protons in a cooler ionized hydrogen plasma (with temper- 
ature indicated along the curves). Solid line: fully relativis- 
tic expression of Stepney and Guilbert (1983), dash-dotted 
line: the non relativistic expression from Spitzer's (1962) 
theory. Dashed line: Spitzer's formula, but averaged over 
a relativistic instead of a classical Maxwellian. Spitzer's 
treatment is accurate to better than 5% for proton tem- 
peratures < 100 MeV. 



that any charge imbalance will quickly lead to the buildup 
of a strong restoring electric field. 

Eq. ([!]) is only valid for non-relativistic conditions, 
whose validity needs to be checked for the high proton 
temperatures in an ADAF. Analytic expressions for rela- 



tivistic temperatures exist only in special cases. 3tepney 
(1983| ) has derived such an expression for the rate of 



Tal Te" I gives some values for the estimated stopping Guilbert (1983 ) 



transfer of energy between populations with relativistic 
Maxwellian distributions in terms of an integral over the 
scattering cross-section. In the case of hot protons heating 
cooler electrons the Rutherford cross-section is the rele- 
vant cross-section. An expression in closed f orm for the 
heating rate (in ergs cm~ 3 sec -1 ) is given by Stepney fc 



depth | from Eq. (|6|)| and the calculated stopping depth 
[according to Eq. ([|)] by following a penetrating proton 
through an isothermal electron layer until it has lost all of 
its energy in excess with respect to the ambient electrons. 
The estimated stopping depth is generally higher than the 
numerical result due to the reasons mentioned above. For 
low electron temperatures these discrepancies are quite 
conspicuous. 

The stopped protons accumulate at their stopping 
depth in the neutron star atmosphere. To maintain charge 
neutrality, an equal number of electrons has to move to the 
same location. The ADAF electrons entering the disk are 
stopped over a much short distance, compared with the 
protons. To make them move to the location where the 
protons come to rest in the disk, an electric field has to 
develop in the neutron star atmosphere, such that the re- 
sulting electron current has just this property. We assume 
that such an equilibrium actually develops, on the grounds 



<lstp 



3m c (fcT c -fcTp) , » 

2m p Uc U P CTT C X K 2 {\/6 B )Ki{l/e p ) mA 



2(9o+9 p ) 2 + 



(7) 



where 9 = kT e /m e c 2 ,# p = kT p /m p c 2 , and Kq, Ki, 
K2 are Bessel functions. 

The proton-electron heating rate according to the 
stopping formula derived from Spitzer's theory is given 

by 



i p v p s(E p ) 



(8) 



By averaging Eq. (J8j) over a thermal velocity distri- 
bution we obtain a rate that can be compared directly 
with Eq. (15) from Stepney & Guilbert. We do this av- 
eraging both with a classical Maxwell distribution and a 
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relativistic Maxwell distribution. The results of this com- 
parison are shown in Fig. [lj. The heating rates are almost 
indistinguishable up to proton temperatures of roughly 
100 MeV. At higher proton temperatures the deviations 
increase steadily. But for the proton temperatures used in 
our simulation (kT p < 100 MeV) Spitzer's formalism is 
accurate to better than 5% compared to Stepney's result. 
Therefore we are confident that the usage of the classical 
stopping formula following [gpitzer (1962 ) together with 
a classical Maxwell distribution for the protons will yield 
the correct results for to energy transfer from virialized 
protons to electrons near a neutron star surface. 

2.2. The model atmosphere 

We compute the density distribution through our plane- 
parallel one dimensional model atmosphere from the equa- 
tion of hydrostatic equilibrium, 



dP 



_9_ 



dP P (r) 
dr 



(9) 



Here g = G ^l* is the gravitational acceleration, which 
we treat as constant as the atmosphere is thin compared 
to the neutron star radius, and K es = 0.40 cm 2 g _1 is the 
electron scattering opacity, ^^-(r) denotes the force ex- 
erted by the deceleration of the protons as a function of 
optical depth. The change of momentum (exerted on the 
atmosphere) of an incoming proton with respect to the 
vertical can be expressed by 



, Av cos ( 



(10) 



The force dP v /dr is evaluated by recording the change 
of velocity Av of the incoming protons as a function of 
optical depth. Since we only study cases with luminosities 
below the Eddington limit we neglect the contribution of 
the radiative pressure. 

We solve our pressure profile by starting with tq=0 at 
the top of our atmosphere and integrate Eq. ( |Io| ) to the 
maximum optical depth of our model atmosphere r max . 
We usually set r max = 10. At this optical depth the energy 
flux of the protons is already negligible. 

2.3. Accretion from the ADAF 

The protons from the ADAF have a temperature of the 
order of the virial temperature [Eq. ([!])]. As mentioned 
above we use a Maxwell distribution for the proton veloc- 
ities. Additionally, the ADAF has a rotation rate, which 
is somewhat smaller than the Keplerian rate fix (e.g. 
Narayan & Yi 1995). When computing the penetration of 
the protons we take the Maxwellian velocity component of 
the protons. Further we take into account the component 
of velocity tangential to the neutron star surface due to 
the ADAF rotation. Instead of using a detailed model of 
an ADAF in which these velocity components can in prin- 
ciple be determined quantitatively, we model it with a 



temperature T = £T V i r and tangential velocity v — rjR^l^, 
i.e. we use the parameters £ and n to scale these quanti- 
ties. Both the tangential and vertical velocity components 
add to the energy deposited, but are not equivalent since 
the thickness of the heated layer is smaller the more the 
protons enter tangentially. 

2.4. The radiative transfer 

The radiative transfer equation appropriate for our model 
is 



(ii) 



where is the intensity as a function of frequency v 
and photon angle fi = cos(#). The coordinate z is mea- 
sured from the top of the atmosphere where the density 
is virtually zero (z — 0) downwards into the neutron star 
atmosphere until z = z max (corresponding to r max ). 

The bremsstrahlung absorption coefficient 0$ is given 

by 



oft = 3.7 x 10 s n c «i 



I _ e -hv/kT 



T 



■9ff 



(12) 



in units of cm . The symbol gg represents the velocity 



averaged Gaunt factor for bremsstrahlung [see e.g. Rybicki 
|fc Lightman (1979 )]. The Compton extinction coefficient 
a% is simply 



n e o„ 



(13) 



where o c v is the electron cross section with Klein-Nishina 
corrections. For the non-relativistic limit this equals the 
Thompson cross section cry. 

The bremsstrahlung emission jff is given by 



ff = 5.4 x 10" 



30 



-hv/kT 



n c n { 



T 



■9ff 



(14) 



in units of erg cm -3 sec -1 ster -1 . The Compton 'emissiv- 
ity' (i.e. radiation scattered to p, v from some other //, v') 
is 



(J^j Ifi'y cr c (fl', v'—>fJ,, v) dp!dv' 



(15) 



Here a c ([jf,v' — * (J-,v) is the Compton Scattering 
Kernel (CSK). An efficient algorit hm for run-time com- 
put ation of the CSK is described by Kershaw et al. (1986 ) 
and Kershaw (1987 ). Both er£ and <r c (/i', v' — v) depend 
on the electron temperature. 

Whenever the temperatures of our solutions are low 
enough, kT m e c 2 , we can simplify the Compton emis- 
sivity Eq. (|l^) without compromising accuracy. We intro- 
duce the angle-averaged Compton cross section a c (y' ' —*v) 
defined as 



(16) 
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A semi-analytic expression for a c (y' — ► v) is given by 
Poutanen & Svcnsson (199E). The angle-averaged inten- 
sity J„ is defined as 



Ju 



(17) 



By replacing <r c (//, z/— >//, v) with its average a c iy' ' —* 
v) in Eq. (|l^) one obtains an approximative expression for 
the Compton emissivity 



(18) 



f v = n j f^j J v ,a c (y' -^v)dv' , 



which is independent of /i, and which involves only a single 
integral instead of the double integral of Eq. (filf). 

Eqs. (0, II, H, ||, constitute the complete set 

of equations for the radiative transfer. Since Eq. (11) and 
Eq. (|lE / 18 ) are mutually dependent, the solution cannot 
be found by direct evaluation. To solve the system, we 
use a standard Lambda Iteration sche me [see e.g. Rutteij 
(1999|)"|and references therein ; see also Zanc ct al. (199(f h 
Poutanen fc Svcnsson (1996| )]. 

We discretize frequency v with equal spacing in log 
and photon angle /i = cos(6>) according to the roots of the 
Legendre polynomials of order n. At the start of the pro- 
cedure we choose j't „ = 0, or any other initial guess that 
might be appropriate. For each Vi and /ij we now integrate 
the transfer equation (Eq. [Tlj ) from z = to z = z m ax for 
Hj > 0, or from z = z max to z = for Hj < 0. At z — we 
impose zero inward flux as boundary condition, while at 
z — Zmax we choose a Planck function at temperature Tbc 
as the starting value for the integration. The temperature 
Tbc of this Planck function is chosen such that: 



oTb C = 2n I fid^ I dv I^z = z max ) 



(19) 



so that at the bottom of the atmosphere the net flux is 
zero. 

After performing this operation for every Vi and fij , we 
can evaluate the new j° by employing Eq. (|l5|/|l8|). This 
new ffy is inserted into Eq. (|ll]) and the integration of 
Eq. ( |ll| ) is repeated along the lines sketched above. After 
each iteration the relative difference of the solution with 
the previous iteration step is computed. If this relative 
difference drops below 10 -3 , we assume that the solution 
has converged. Since in our solutions the scattering optical 
depths remains always of the order of a few, we need not 
worry about the well-known convergence problems of the 
lambda iteration procedure sketched above. 

Once the radiative transfer solution is obtained, the 
radiative heating and cooling rates can be evaluated by 
flux differences. The flux at each position x is defined as 



F(z) = 2tt / /idfi dv I^^ix) 



The net cooling rate is then 



dF(z) 
dz 



(20) 



(21) 



The radiative transfer is solved on a grid set up in 
terms of Thomson optical depth when solving the pres- 
sure profile (Sect. [^). This grid is then transferred to 
the vertical z coordinates for the radiative transfer code 
by the usual 



z(r) 



(n c (r)CT T ) l dr. 



(22) 



Thus we immediately obtain the net cooling rates 
A^ d (r) as a function of Thomson optical depth. 

We also checked our solution for the importance of pair 
processes, calculating the pair number densities according 
to Zanc et al. (1998). We find that in our model the elec- 



tron temperatures are too low and the electron densities 
too high for pair production to become important. 

2.5. Energy balance from heating and cooling 

We start our calculations with an isothermal atmosphere 
in hydrostatic equilibrium according to Sect. |2.2| . For 
the initial electron temperature we set T e = 1 keV 
throughout the layer. After each time step we obtain the 
heating rates A+(r) from the Coulomb interactions and 
the cooling rates A^ d (r) due to the radiative processes 
bremsstrahlung and multiple Compton scattering as a 
function of optical depth. The time step of every cycle in 
our simulation is adjusted to the shortest energy exchange 
time scale occurring in the calculation. 

Additionally we include the energy redistribution due 
to electron thermal conductivity. Generally the flow of 
heat per unit area, Q, in the presence of a temperature 
gradient, is given by (see e.g. Bpitzer 1962 ) 

Qw-kVT. (23) 



by 



The heat conductivity k for a fully ionized gas is given 



1.85 x 10" 



rp5/2 

TnT 



erg sec 



(24) 



Thus the change of enthalpy dw/dt due to electron 
conductivity can be expressed by the divergence of the 
heat flux Q, 



dw dT 

^=pc p -=V.(.VT) 



A r 



(25) 



where c p is the specific heat at constant pressure. 

Now we can calculate the total change of enthalpy per 
time step, 



Aw(t) AT(t) 
- pc p - 



A ?M+ A ra»+ A cond(T). (26) 



At r " At 
With Eq. (^6|) we are able to calculate the change of 
temperature and thus the new temperature profile after 
each time-step. With the new temperature profile we can 
update the hydrostatic structure according to Sect. |2.2[ 
We follow the simulation until the Coulomb heating is bal- 
anced by the radiative cooling and energy redistribution 
by conductivity and the temperature and density struc- 
ture as well as the spectra have reached a stationary state. 
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Fig. 2. upper left: emergent spectrum 
(solid line), thermalized downward di- 
rected spectrum (dashed line) and 
blackbody input spectrum at the base 
(dotted - dashed line); upper right 
panel: combined heating rates from pro- 
ton heating and electron thermal con- 
ductivity (solid line) and radiative cool- 
ing rates due to Comptonization and 
bremsstrahlung (squares) - the dotted 
line shows the rates from electron con- 
ductivity alone; lower panels show from 
left to right and from top to bottom 
the electron temperature T a , pressure 
P, electron density n e and the geomet- 
rical depth of the layer z, for a solution 
with F v = 1.6 x 10 24 erg cm -2 sec -1 , 



optical depth t 



1 .0 

optical depth r 



±T vi r and O. 



3. Results of the model computations 

For our model calculations we use fixed values for the mass 
M* = IMq and the radius i?* = 10 6 cm for the neutron 
star. For the rotational angular velocity of the accreting 
material we use rj = 0.5. The results depend on the pro- 
ton energy flux per unit area of the neutron star surface. 
In order to make the results interpretable in terms of a 
total luminosity, we assume that 60% of the neutron star 
surface is involved in the accretion process. This takes 
i nto account the approxim ate vertical extent of an ADAF 
(|Popham fc Sunyacv 200ll) . 



cooler than the "hot" solutions from Iurolla et al. (1994) 



3.1. Dependence on proton energy flux and 
temperature 

Fig. ^ shows an equilibrium solution for a proton energy 
flux F p = 1.6 x 10 24 erg cm -2 sec -1 ,^ = rj =0.5, yielding a 
luminosity of = O.lLEdd- The neutron star atmosphere 
is clearly divided into two parts. A hot surface layer is 
separated by a sharp temperature front from the much 
cooler bottom part. The maximum temperature in the hot 
part of the atmosphere is T c « 50keV. This is consider- 



ably hotter than previous calculatio ns showed (e.g. |Almc 
fc Wilson 1973) ; |Turolla et al. 1994|) but also significantly 



and Zane et al. (1995). We did not find any comparable 



hot solutions. 

Fig. H also shows that energy redistribution through 
electron conductivity gets important at the transition zone 
from the hot to the cool part. The emergent model spec- 
trum in Fig. [2] illustrates that the hot part acts as an ef- 
fective Comptonization layer with an optical thickness of 
order unity. The downward directed photons are almost 
completely thermal ized in the cold part. This justifies our 
assumption of Sect. 2A that the downward flux of photons 
is converted into a black body flux of upward soft photons 
of temperature Tbc at the base of our model since almost 
all the thcrmalization has already taken place. 

Fig. U shows the emergent spectra and temperature 
profiles through the atmosphere at a fixed proton energy 



flux level F p = 1.6 x 



10 24 erg cm 2 for various proton 



temperatures, £XUr- The higher the initial proton tem- 
perature the lower is the electron temperature in the hot 
part of the atmosphere and vice versa for the cold part. 
The optical depth of the hot part is of order unity in every 
case. As a consequence harder spectra are produced with 
lower initial proton temperatures. This can be explained 
with the decreasing penetration depth of the protons at 
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Fig. 3. Emergent model spectra and temperature profiles at a fixed proton energy flux F p = 1.6 x 10 24 erg cm -2 sec -1 
for various initial proton temperatures, £T v i r , with £=0.25,0.5,0.75, and ?7=0.5. Higher proton temperatures produce 
lower electron temperatures in the hot part of the atmosphere and softer emergent spectra. The dotted dashed line 
shows the corresponding blackbody spectrum with the effective temperature of the layer, T q = (Fp/er) 1 / 4 . 




Fig. 4. Emergent spectra and temperature profiles for various values of the proton energy flux (in erg cm -2 sec -1 , 
denoted by the numbers of the lines) at constant proton temperature, £ = 0.5. With decreasing energy flux the optical 
depth of the heated layer shrinks. 



smaller velocities: the heating rates in the upper part of 
the atmosphere are increased and thus the temperature is 
higher there. 



trum with respect to the blackbody of the layer effective 
temperature at the lowest luminosity. 



Fig. |^ and Fig. |5| show the dependence of the solution 
on the proton energy flux at constant proton temperature, 
£ = 0.5. At low proton energy flux (Fig. ||) the optical 
depth of the hot part is small (r^ot < 10~ 2 ). The solutions 



for these cases are comparable to those found by Zampicri 
et al. ( |1995| ). But the temperature of the hot part in our 
solutions turns out to be much higher than in their solu- 
tions. This hot part causes an (energetically unimportant) 
bump in the emergent spectra at high frequencies due to 
op tically thin brcmsstra hlung emission in this layer. As 
in Zampicri et al. (1995 ) the emergent spectra are signif- 
icantly harder than the blackbody at the layer effective 
temperature, T e g = (Fp/cr) 1 / 4 . We also see a harder spec- 



By further increasing the proton energy flux the tem- 
perature jump moves inward to greater optical depth 
reaching order unity at the highest flux levels. Thereby 
the temperature in the hot part is roughly kept at a con- 
stant value in contrast to the temperature of the cold part, 
which increases continously. At F p = 1.6 x 10 24 erg cm -2 
sec -1 the hot part acts as an effective Comptonization 
layer for the cool thermal photons from the interior. In 
addition the electron thermal conduction smoothes the 
temperature profile and the shock is not as sharp as in 
the low energy flux cases. 
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Fig. 5. Emergent spectra and temperature profiles at low proton energy flux (in erg cm -2 sec -1 , denoted by the 
numbers of the lines) and constant proton temperature, £ = 0.5. The hot part is extremely thin in terms of Thomson 
optical depths. The dotted dashed lines show the corresponding blackbody spectra with the effective temperature of 
the layer. 



3.2. A reverse photosphere effect 

Evident in our model spectra is an excess of soft pho- 
tons (see e.g. Figs. Hflll) w hh respect to a blackbody 
with the effective temperature of the corresponding energy 
flux. Though energetically unimportant it is an interest- 
ing feature of our model spectra. This soft photon excess 
is caused by the reverse temperature profile of the atmo- 
sphere, i.e. the high temperatures are above the cooler 
parts (in contrast to e.g. the solar photosphere). Fig. || 
shows the absorption optical depth as a function of the 
photon energy and electron temperature for an equilib- 
rium solution with £ = 0.5 and r\ = 0.5. The line for 
r a = 1 is emphasized. For low photon energies r a = 1 
is found high in the neutron star atmosphere, where the 
temperature is high. For high frequency photons the atmo- 
sphere gets more and more transparent and T a = 1 occurs 
at lower temperatures. Thus the high frequency photons 
decouple at lower temperatures and propagate outwards. 
As a result there is an excess of low energy photons, which 
decouple from the atmosphere at higher temperatures. 

The same arguments explain the harder spectra at 
low proton energy fluxes. Here the temperature profile in 
the atmosphere is dominated by the increase of temper- 
at urc with optical depth , as has already been explained 
by Zampicri ct al. (1995 ). Thereby the high energy pho- 
tons decouple in the deeper, hotter layers and propagate 
out freely This contribution to the high energy part of 
the spectrum causes the harder spectrum compared to the 
blackbody of the layer effective temperature. 



4. Discussion and conclusions 

The main conclusion from our work is that a neutron star 
surface embedded in a hot ADAF-type accretion flow does 
not act like a simple blackbody thermalizer. Instead the 




0.00 



0.100 1.000 
photon energy [keV] 



00.000 



Fig. 6. Absorption optical depth as a function of photon 
energy and atmospheric temperature for an equilibrium 
solution with £ = 0.5 and r\ — 0.5. The line depicting 
r a = 1 is emphasized. 

interaction with the hot ions produces spectra with pro- 
nounced high energy, Comptonized tails. 

We have considered the accretion of protons onto the 
surface of a neutron star for a wide range of proton energy 
fluxes. In our model the stopping depth of the protons 
is obtained self consistently. We allow for the scaling of 
the proton temperature T p = £T v ir, which in our model 
controls the vertical velocity component of the protons. 
Further we consider accretion from an ADAF, where the 
protons have a tangential velocity component v = r/Q^- 
By variation of the proton temperature and the proton 
energy flux a range of X-ray spectra can be produced. Our 
model computations show that the local proton energy 
flux determines the optical depth of the hot part whereas 
the temperatures through the atmosphere are determined 
by the initial proton temperatures. 
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Previous models considered spherical accretion onto 
the neutron star. The free infall velocity of the protons 
was used in such models. If accretion takes place from 
an accretion disk or an ADAF, the vertical velocity com- 
ponent is considerably smaller than than the free infall 
velocity. This has important consequences on the emer- 
gent spectra as our solutions suggest: at the same local 
energy flux a decrease in the vertical proton velocity in- 
creases the temperature of the hot part of the atmosphere, 
and X-ray spectra similar to those of neutron stars in their 
hard states can be produced. Fig. [7] shows a comparison of 
a model spectrum with observed spe ctra of X-ray burster s 
in their low/hard states (taken from |Gilfanov et al. 1998) ). 

A question left unanswered by our studies is how the 
spectrum from the surface is modified by the overlying 
accretion flow. E.g. an overlying optically thick boundary 
layer might significantly thermalize hard photons from the 
neutron star surface again. But we can conclude that the 
input spectrum at the bottom of a boundary layer is not 
Planckian at high proton energy fluxes. If the boundary 
layer is optically thin the contribution to the observed 
hard spectral component at high proton energy fluxes 
might be quite important. A discussion whether there is an 
optical ly thi n boundary layer can be found in e.g. Barr 
et al. ( 2000 ). For low accret ion rates the bounda ry layer 
is indeed optically thin (e.g. King fc Lasota 1987 ). 
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Detailed studies about the accretion disk boundary 



layer around a NS have been performed by Popham & 



Sunyacv (2001). They show that the accretion disk near 
the NS is radially and vertically extended and that the an- 
gular velocity f2 is reduced with respect to the Keplerian 
value. Such a boundary layer naturally produces a pro- 
ton illumination scenario as described in our model. Their 
boundary layer is optically thin to absorption. For high 
luminosities this region has a radial extent larger than 
one stellar radius, i.e. the hard radiation from the NS sur- 
face, produced by proton illumination, can easily propa- 
gate outwards without being thermalized. The boundary 
layer even further enhances Comptonization as the gas 
there is hot (> 10 8 K). 

Our results emphasize the importance of the neutron 
star surface for the contribution to the hard spectral com- 
ponent in low mass X-ray binaries (LMXBs). A prereq- 
uisite for the occurrence of a hard spectral component 
is a high local proton energy flux. The local proton en- 
ergy flux onto the neutron star surface depends on both 
the mass accretion rate and the size of the accretion belt 
on the neutron star. Our model computations have local 
validity, i.e. if the involved accretion belt is small, high 
proton energy fluxes on the neutron star surface can even 
occur at accretion rates well beneath the Eddington limit. 
Therefore our model does not exclude the possibility of 
hard spectral components from the neutron star surface 
at low luminosities, as is observed. 
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Fig. 7. Comparison of an emergent model spectrum with 
£ = 0.25 and rj — 0.5 with observed neutron star spectra 
in their hard states. The spectra are arbitrarily shifted in 
the vertical direction for clarity. 
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